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Abstract. We present continuum observations of a small 
sample of high-redshift, radio- loud AGN (radio galaxies 
and quasars) aimed at the detection of thermal emis- 
sion from dust. Seven AGN were observed with IRAM 
and SEST at 1.25 mm; two of them, the radio galax- 
ies 1243+036 (z ~ 3.6) and MG 1019+0535 (z ~ 2.8) 
were also observed at 0.8 mm with the JCMT submil- 
limetre telescope. Additional VLA observations were ob- 
tained in order to derive the spectral shape of the syn- 
chrotron radiation of MG 1019+0535 at high radio fre- 
quencies. MG 1019+0535 and TX 0211-122 were expected 
to contain a large amount of dust based on their depleted 
Lya emission. The observations suggest a clear 1.25-mm 
flux density excess over the synchrotron radiation spec- 
trum of MG 1019+0535, suggesting the presence of ther- 
mal emission from dust in this radio galaxy, whereas the 
observations of TX 0211 — 122 were not sensitive enough 
to meaningfully constrain its dust content. On the other 
hand, our observations of 1243+036 provide a stringent 
upper limit on the total dust mass of < 10 s Mq. Finally, 
we find that the spectra of the radio-loud quasars in our 
sample (2 < z < 4.5) steepen between rest-frame radio 
and the far-infrared. We discuss the main implications 
of our results, concentrating on the dusty radio galaxy, 
MG 1019+0535. 

Key words: Galaxies: active - Galaxies: fundamental 
parameters - Galaxies: ISM - ISM: dust - Radio contin- 
uum: galaxies 
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1. Introduction 

In recent years several attempts have been made to study 
the interstellar medium in high-redshift galaxies. Promis- 
ing results have been obtained by millimetre and submil- 
limetre observations, which sample the rest-frame far-IR 
emission of galaxies at z > 2 (Andreani et al. 1993). The 
study of dust in distant galaxies is very important be- 
cause it provides information on the physical state of the 
ISM, and its relation to other properties of the galaxies, 
such as activity and evolution. For example, if it is cor- 
rect to assume that the rest-frame far-IR luminosity Lfir, 
is a measure of the star-formation rate, then this could 
provide a way to estimate the evolutionary state of high-z 
galaxies. Establishing whether dusty galaxies are common 
at high redshifts is also important in order to evaluate the 
effects of dust obscuration on surveys of quasars and pro- 
togalaxies carried out in optical bands (Smail, Ivison & 
Blain 1997). 

Several active galaxies with z > 2 have been de- 
tected in submillimetre/millimetre bands, suggesting the 
presence of large amounts of dust in their host galaxies 
(M d - 10 8 " 9 M Q ; see Hughes, Dunlop & Rawlings 1997 
for a recent review; Omont et al. 1996b; Ivison et al. 1998). 
At high redshift, dust is also thought to be present in 
damped Lya absorption systems (Pettini et al. 1994; Pei, 
Fall & Bechtold 1991), in very red galaxies (Hu & Ridg- 
way 1994), and in high-z radio galaxies (Cimatti 1996 and 
references therein). A substantial amount of dust is also 
expected in theoretical models of the evolution of galaxies 
at high-z (Mazzei & De Zotti 1996 and references therein). 
Finally, it is important to recall that molecular gas has 
been observed in a few distant active galaxies, allowing 
a direct estimate of the dust-to-gas mass ratio at large 
cosmological distances for the first time (see Omont et al. 
1996a; Ohta et al. 1996). 

We recall that particular caution is needed in the inter- 
pretation of the submillimetre/millimetre data of radio- 
loud objects. The existence of a thermal dust emission 
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excess should be confirmed by checking whether the sub- 
millimctrc/millimctrc flux densities represent a real ex- 
cess over the synchrotron spectrum. The radio galaxy 
B2 0902+343 (z ~ 3.4) is an example of the problem, 
where the observed 1.3- mm flux density is not due to 
thermal emission from dust, but to the tail of the radio 
synchrotron spectrum (Downes et al. 1996 and references 
therein). Therefore, high-frequency (~ 8 — 43 GHz) radio 
observations are extremely important to derive the shape 
of the synchrotron spectrum and to estimate its contribu- 
tion to the millimetre region. 

In the present paper we present the results of milli- 
metric observations of a small sample of radio galaxies 
and radio- loud quasars with z > 2, and additional high- 
frequency VLA observations of MG 1019+0535, a radio 
galaxy at z ~ 2.8. Throughout this paper we assume 
H = 50kms _1 Mpc -1 , q = 0.5 and define /150 = H o /50. 

2. Observations and data reduction 

We have observed seven objects with 2 < z < 4.5; four ra- 
dio galaxies and three radio-loud quasars. Table 1 shows 
the target details. The radio galaxies TX 0211 — 122 and 
MG 1019+0535 were observed because their UV spectra 
are suggestive of extinction by dust (van Ojik et al. 1994; 
Dey et al. 1995), MRC 0943-242 because of the the pres- 
ence of a large amount of neutral hydrogen surrounding 
the host galaxy (Rottgering et al. 1995) and 1243+036 
because of its very high redshift (van Ojik et al. 1996). 
Among the quasars, we recall that PKS 1251— 407 is, to 
date, the most distant known radio-loud quasar (Shaver, 
Wall & Kellermann 1996). MRC 1043-291 (Kapahi ct al. 
1997) and PKS 1354-107 (Shaver et al. 1997, in prepa- 
ration), were selected because of their convenient observ- 
ability during the observing run. 

2.1. SEST observations 

The observations with the 15-m SEST (Swedish— ESO 
Submillimetre Telescope) were made at La Silla, Chile, 
during 1995 July and 1996 May. The telescope was 
equipped with a 3 He-cooled 1.3- mm bolometer, with a 
central frequency of 250 GHz and a bandwidth of around 
50 GHz. The bcamsizc is about 24" (FWHM) and the 
typical sensitivity is 200mJys~ 1 / 2 . The observations were 
carried out during nights with rather low opacities. The 
opacity was checked by several skydips taken during the 
nights (i.e. the telescope was moved to six different el- 
evations where the bolometer integrated for 10 s and a 
calculation of the zenith opacity was made). The pointing 
was checked using quasars as astrometric reference sources 
and the absolute flux calibration was achieved by observ- 
ing the planet Uranus. We estimate that the flux calibra- 
tion has a typical uncertainty of 20%. The observations 
were made in beam-switching mode with a beam throw of 
70", and the data were reduced according to method de- 



scribed by Andreani (1994) (see also Andreani et al. 1993). 
We have found a linear decrease of the standard deviation 
with t~ x / 2 , indicating that the observing conditions were 
stable and limited by the sky noise. 

2.2. IRAM observations 

Observations with the IRAM 30-m antenna were made in 
1996 March with the MPIfR 1.3-mm 7-channel bolometer 
array. The bolometers are separated by 22" in an hexag- 
onal arrangement surrounding the central pixel. The in- 
dividual beamsize is about 11" (FWHM), and the typ- 
ical sensitivity of eOmJys" 1 / 2 . The ON-OFF technique 
uses the wobbler with a beamthrow of 33", integrating for 
10 s per subscan. The opacity was monitored every 1 — 
1.5 hr with a skydip procedure, and was around 0.4 during 
the first night and 0.2 during the second. In both cases, 
the atmosphere was quite stable. The flux calibration was 
achieved by observing Uranus, and the uncertainties are 
again of order 20%. The fluxes of the six outer channels 
were averaged and subtracted from the central channel. 

2.3. J GMT observations 

The radio galaxies 1243+036 and MG 1019+0535 were 
also observed at the JCMT (the 15-m James Clerk 
Maxwell Telescope, Mauna Kea, Hawaii) using the single- 
element 3 He-cooled UKT14 bolometer (Duncan et al. 
1990), coupled with a broad-band filter (v = 384 GHz; 
Av = 30 GHz FWHM). The observations were made us- 
ing a 65-mm focal plane aperture, resulting in a FWHM 
beamwidth of 16.5". Sky emission was subtracted by chop- 
ping the secondary mirror in azimuth, at a frequency of 
7.8 Hz, with a throw of 60". 

The radio galaxy 1243+036 was observed on 1996 
April 26 during excellent weather conditions. Calibra- 
tion was taken from Mars (which set near the beginning 
of the shift) and the secondary calibrators NGC 2071IR, 
CRL 618, IRC +10216 (CW Leo) and 16293-2422. 3C 273 
was used as a pointing source and bootstrap calibra- 
tor (observed eight times). The 384-GHz zenith optical 
depth was typically around 0.29, but it appeared to de- 
crease towards the end of the shift (0.23). In total, 880 
pairs of 16 s (i.e. 8 s in each beam) were obtained. The 
raw, uncalibrated dataset gave a S/N of 0.99. After de- 
spiking, calibration and statistical testing, the final re- 
sult was 7.0 ± 3.1 mJy The calibrated 800-/im flux of 
3C273 was 8.2 ± 0.1 Jy; IRC +10216 (which is variable) 
was 5.4 + 0.1 Jy. 

The radio galaxy MG 1019+0535 was observed over 
a period of 5.5 hr during 1996 April 24. A total of 3hr 
was spent on-source, the remaining 2.5 hr being devoted 
to calibration. Flux calibration was provided by regular 
measurements of IRC +10216, as was a determination of 
the 375-GHz zenith opacity. Observations of local pointing 
sources kept rms pointing errors below 3" and associated 
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Table 1. Summary of the observations. 



Name 


Class 


z 


Tel 


Date 


^obs 


Arest 


tint 


T 


S v 


S v /a 












(H 


H 


(s) 




(mJy) 




TX 0211-122 


RG 


2.34 


SEST 


13 July 1995 


1270 


380 


4500 


0.19-0.22 


-0.53+3.99 


- 


MRC 0943-242 


RG 


2.93 


SEST 


06 May 1996 


1270 


323 


5500 


0.11-0.13 


4.33+3.29 


1.3 


MG 1019+0535 


RG 


2.76 


IRAM 


19-20 March 1996 


1250 


332 


5760 


0.30-0.40 


2.13+0.47 


4.5 


MG 1019+0535 


RG 


2.76 


JCMT 


24 April 1996 


790 


210 


10800 


0.30-0.40 


14.70+4.60 


3.2 


MRC 1043-291 


RQ 


2.13 


SEST 


07 May 1996 


1270 


406 


800 


0.15-0.17 


29.38+8.15 


3.6 


1243+036 


RG 


3.58 


JCMT 


25 April 1996 


790 


172 


14080 


0.23-0.29 


7.00+3.10 


2.3 


1243+036 


RG 


3.58 


IRAM 


19-20 March 1996 


1250 


273 


2040 


0.20-0.40 


2.61+0.86 


3.0 


PKS 1251-407 


RQ 


4.46 


SEST 


07 May 1996 


1270 


233 


3800 


0.14-0.19 


8.00+3.1 


2.6 


PKS 1354-107 


RQ 


3.00 


SEST 


06 May 1996 


1270 


317 


2100 


0.14-0.16 


10.92+5.11 


2.1 



Note: RG — radio galaxy; RQ — radio-loud quasar; tint — integration time on the source; r — opacity during the observations. 

flux losses below 10%; the overall uncertainty in the flux 
measurement was therefore dominated by the poor S/N. 
After careful editing, the MG 1019+0535 measurements 
indicated a 375-GHz flux density of 14.7 ± 4.6 mJy — a 
marginal detection. 

2.4- VLA observations 

MG 1019+0535 was also observed at the Very Large Array 
(VLA), New Mexico. Radio continuum observations in the 
X (8.3 GHz), U (15.0 GHz) and K (23.0GHz) bands were 
carried out in the D configuration on 1996 July 19 (with 
a bandwidth of 50 MHz) , and in the Q band using the A 
configuration on 1996 December 26. 

The total on-source time was 10 minutes each in the X 
and U bands, and 20 m in the K band. The visibility aver- 
aging time was 10 s in all bands. The phase was calibrated 
with calibrator source 1024—008 which was observed every 
10 m, and 3C286 was observed as a flux calibrator. The 
data were reduced using standard AIPS procedures. The 
used synthesized beam sizes, are listed in Table 2. Maps 
with a cell size of 0.5" were produced for each band. Maps 
were made and CLEANed using the MX routine within AIPS. 
The radio source was unresolved — no evidence was found 
for any extended emission. The total flux density, as mea- 
sured from the maps in each band, is also given in Ta- 
ble 2. The major source of uncertainty for the measured 
flux density is the uncertainty in the flux scale, which we 
estimate to be of order 10%. 

The Q band observation was carried out over a period 
of 1.7 hr using twelve VLA antennas. The total bandwidth 
for the observations was 100 MHz, centred at 43.34 GHz 
(6.9mm). Again, the observing and calibration proce- 
dures were standard. After checking the pointing accu- 
racy of the antennas, brief observations of MG 1019+0535 
were sandwiched and interspersed with measurements of 
1055+018, a bright, compact calibrator. The flux density 



of the galaxy was tied to that of the calibrator (4.62 Jy) 
which, in turn, was tied to the flux density of 3C 286 
(1.49 Jy on the shortest baselines). The noise level at the 
position of the galaxy agreed well with that measured us- 
ing the calibrated visibilities, giving a 3cr upper limit of 
1.94 mJy. 



Table 2. VLA observations of MG 1019+0535. 



Date 


fobs 


frest 


Flux 


Config 


Beam size FWHM 


(1996) 


(GHz) 


(GHz) 


(mJy) 


& Band 


(") 


19 July 


8.3 


31.2 


52.7+5.3 


D,X 


7.36x5.97 


19 July 


15.0 


56.4 


23.2+2.3 


D,U 


4.31x3.64 


19 July 


23.0 


86.8 


9.3+0.9 


D,K 


2.80x2.52 


26 Dec 


43.3 


162.8 


3a <1.94 


A,Q 





3. Results 

Table 1 shows the flux densities measured for our targets. 
The radio galaxy MG 1019+0535 is the only source de- 
tected at submillimetre wavelengths at greater than the 
4-cr significance level. In this section we discuss the main 
implications of our observations for each individual source. 
The dust masses are estimated using the formula: 

M d = ? - * , (1) 

where S is the flux density, and z/ rest are, re- 

spectively, the observed and rest-frame frequencies, 
£>l is the luminosity distance, B is the black-body 
Planck function, T<j is the dust temperature and «d = 
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0.67(^ res t/250 GHz) 2 cm 2 g _1 is the adopted mass absorp- 
tion coefficient. In order to obtain limits on Md, it is nec- 
essary to assume a temperature for the grains. We as- 
sume a representative temperature Td=60K, consistent 
with the typical temperatures estimated for high-redshift 
active galaxies. For a discussion about the uncertainties of 
dust masses and temperatures, see Hughes et al. (1997). It 
should be emphasised here that, besides the uncertainties 
on Xd and Kd, the absolute dust masses are also strongly 
dependent on the choice of Hq and qo; in fact, they change 
by a factor of 4 for values of H ranging from 50 to 
100 kms -1 Mpc~ x , and by a factor of around 2 if qo is 
changed from go = 0.1 to qo = 0.5. n 

For each source we have searched NEEu in order to de- 
rive the spectral energy distributions (SEDs) over a broad 
range of frequencies. If the error of the flux density is un- 
known, we assume an uncertainty of 5%. Whenever the 
significance of the flux density is < 3c, we provide an up- 
per limit at the 3-cr level. Whenever more than one flux 
density value is available at the same frequency, we plot 
all the available values. Figs. 1 and 2 and Table 3 show 
the SEDs derived for our targets. 

3.1. Notes on individual sources 
TX 021 1-122 

The optical emission line ratios of this radio galaxy are 
highly anomalous, with the flux of the Lya line relatively 
weak, and that of the N VA1240 line relatively strong com- 
pared to those of CivA1549, HenA1640 and Cm]A1909 
(van Ojik et al. 1994). This was interpreted as being due to 
a strong starburst and a large amount of dust. According 
to the 8.2-GHz flux density and spectral index (Rottgering 
et al. 1994; Carilli et al. 1997), the expected synchrotron 
flux density at 1.3 mm is only ^0.3 mJy. Our SEST obser- 
vation yields a 3-cr upper limit of ^12mJy which limits 
the total amount of dust, according to equation (1), to 
< 1.1 x 10 9 Mq. Although this limit is not very stringent, it 
does suggest that the total amount of dust cannot be much 
larger than that of the most dusty known active galaxies: 
for instance, the quasar BRI 1202-0725 has M d ~ 10 9 M 
(Isaak et al. 1994; Hughes et al. 1997). Finally, we recall 
that the amount of molecular hydrogen estimated from 
observations of CO is < lO n M0; not much greater than 
that of nearby gas-rich starburst galaxies (van Ojik et al. 
1997). 

MRC0943-242 

The detection of a halo of neutral hydrogen linked with 
the host galaxy of MRC 0943-242 and the dust which 
might be associated with the neutral ISM (Rottgering et 
al. 1995) prompted the SEST observation of this source. 
We recall that the limit on the molecular hydrogen mass 
is < 1O U M (van Ojik 1995). According to the 8.2-GHz 

1 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, Caltech, under contract 
with the National Aeronautics and Space Administration. 



flux density and spectral index (Rottgering et al. 1994; 
Carilli ct al. 1997), the expected synchrotron flux density 
at 1.3 mm is only ^0.2 mJy. Our 1.3-mm 3-cr upper limit 
of 9.9mJy limits the mass of dust to < 6.8 x 10 8 M Q ). 
Because of the large uncertainty in the total gas content 
(H+H 2 ) (Rottgering et al. 1995; van Ojik et al. 1997), it is 
not possible to meaningfully constrain the dust/gas mass 
ratio in this galaxy. 
MG 1019+0535 

This radio galaxy has spectroscopic properties similar 
to TX 0211 — 122, again indicating the possible presence 
of dust (Dey et al. 1995). Our IRAM observations pro- 
vided a suggestive detection at 1.25 mm, and our JCMT 
observations yielded a marginal detection at 800 /jm. The 
IRAM observations were split into two nights but in this 
case, unlike 1243+036, MG 1019+0535 gave consistently 
positive signal. Although the result is formally significant, 
we consider that our 1.25-mm data provide only a tenta- 
tive detection because of the very weak flux density. We 
note, however, that previous IRAM detections at around 
this level have since proved to be trustworthy — that of 
8C 1435+635, for example (Ivison 1995; Ivison et al. 1998). 

It is important to stress that there are major uncer- 
tainties in the interpretation of the millimetric observa- 
tions of this galaxy. Optical imaging shows the presence 
of two objects separated by about 1.5" — object A, identi- 
fied as the counterpart of the radio source at z = 2.76, and 
object B (Dey et al. 1995). The nature of B is unclear: it 
may be physically related to A, or be a foreground galaxy 
at z ~ 0.66. The problem is that the beam widths of our 
1.25-mm and 800-fim observations include both objects. 
However, if the two objects are unrelated, the depression 
of the Lya line favours component A being the dusty ob- 
ject and the source of the observed flux density at 1.3 mm. 

In order to better constrain the SED of this galaxy, we 
used data from IRAS. However, since MG 1019+0535 is 
not detected, upper limits have been estimated at 12, 25, 
60 and 100 by searching a 1 square degree field centred 
on MG 1019+0535 for sources from the IRAS Faint Source 
Catalogue, adopting the faintest in each band as the up- 
per limit (0.11, 0.17, 0.20 and 0.47 Jy, respectively). This 
crude method relies on the fact that if the FSCs sophisti- 
cated search routines cannot find a point source, then the 
source must be below the 3-cr threshold. The method is 
less prone than some to providing misleadingly low limits 
(Ivison 1995). The implication of this result is dicussed in 
detail in Section 4. 
MRC1043-291 

This radio-loud quasar has radio flux densities of 1.09 
and 0.68 Jy at 408 MHz and 5 GHz, respectively (Kapahi 
et al. 1997). Therefore, if we adopt a spectral index a = 
—0.19 (defined as S v cx v a ), we derive an expected 1.3-mm 
synchrotron flux density of around 326 mJy. Our SEST 
observation provides a flux density around an order of 
magnitude lower than expected, suggesting that the radio 
spectrum steepens rapidly at high frequencies. 
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Table 3. Flux densities. 



Object 



fobs 

(GHz) 


Flux 
(mJy) 




Ref 


Object 


fobs 

(GHz) 


Flux 
(mJy) 




Ref 


0.408 


1020±40 




L81 


MRC 1043-291 


0.408 


1090±60 




L81 


1.465 


189±9.5 




C97 




2.70 


530±27 




WO90 


4.70 


53.8±5.0 




C97 




4.85 


714±51 




W96 


8.20 


24.3±2.5 




C97 




5.00 


320±16 




WO90 


230 


<12 




here 




230 


<24.3 




here 


4.61xl0 5 


2.36+0.24x10" 


-3 


v094 


1243+036 


0.178 


3600±180 




WO90 


0.408 


1050±30 




L81 




0.408 


1810±80 




L81 


1.50 


245±13 




C97 




1.40 


256±13.0 




WB92 


4.70 


56±3 




C97 




2.70 


120±6.0 




WO90 


4.85 


44±11 




Gri94 




4.70 


69.6±7.0 




v096 


8.20 


21±1 




C97 




8.30 


28.2±3.0 




v096 


230 


<9.9 




here 




230 


<2.6 




here 


4.61xl0 5 


5.93+0.59x10" 


-3 


C97 




384 


<9.3 




here 


0.178 


2100±105 




WO90 




1.36x10 s 


1.28+0.13x10" 


-2 


v096 


0.365 


925±31 




D95 




4.61x10 s 


2.84+0.28x10" 


-3 


v096 


0.408 


920±50 




L81 


PKS 1251-407 


0.33 


<400 




S96 


1.400 


454±23 




WB92 




1.40 


260±10 




S96 


1.490 


360±6 




D95 




2.70 


250±10 




S96 


2.700 


180±9 




WO90 




4.85 


238±15 




W94 


4.850 


100±12 




Gri95 




4.85 


237±16 




Gre94 


4.850 


115±17 




B91 




4.90 


200±10 




S96 


4.850 


132±19 




Gre91 




5.00 


220±10 




WO90 


8.333 


52.7±5.3 




here 




15.0 


110+10 




S96 


8.439 


59.1±1.5 




D95 




8.40 


150±10 




S96 


15.00 


23.2±2.3 




here 




230 


<9.3 




here 


23.08 


9.30±0.9 




here 




3.75 x10 s 


2.47+0.25x10" 


-2 


S96 


43.30 


<1.94 




here 




6.82 x10 s 


1.40+0.14x10" 


-3 


v096 


230.0 


2.13±0.47 




here 


PKS 1354-107 


2.70 


260±13 




WO90 


384.0 


14.7±4.6 




here 




4.86 


284±14 




Gri94 


3000 


<460 




here 




5.00 


270±14 




WO90 


5000 


<200 




here 




8.40 


280±14 




S97 


1.20xl0 4 


<170 




here 




8.40 


191±10 




S97 


2.50xl0 4 


<110 




here 




230 


<15.3 




here 


1.36x10 s 


6.40±1.30xl0~ 


■3 


D95 












3.75 x10 s 


1.70±0.17xl0" 


-3 


D95 












4.61x10 s 


7.84+0.40x10" 


-4 


D95 












5.00 x10 s 


5.05+0.50x10" 


-4 


D95 












6.00 x10 s 


4.03+0.40x10" 


-4 


D95 













TX 0211-122 



MRC 0943-242 



MG 1019+0535 



Note: L81 — Large et al. (1981); C97 — Carilli et al. (1997); v094 — van Ojik et al. (1994); Gri94 — Griffith et al. (1994); 
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1243+036 (= 4C 03.24) 

This is the radio galaxy with the highest redshift in 
our observed sample. One spectacular feature is the pres- 
ence of a Lya halo (with a luminosity of ~ 10 44 ergs _1 ) 
which extends over 20" (~ 135 kpc) (van Ojik et al. 
1996). The Lya images, coupled with high-resolution spec- 
tra, indicate that the radio jet is interacting vigorously 
with the gas in the inner region. Perhaps most surpris- 
ing is the low-surface-brightness outer region of the Lya 
halo. Deep spectroscopy shows that it is relatively quies- 
cent (~ 250 km s -1 FWHM), but that there is a veloc- 
ity gradient of 450 km s -1 over the extent of the emission 
(~ 135 kpc). Because the halo extends beyond the radio 
source, it is probable that its kinematics must predate the 
radio source. The ordered motion may be a large-scale ro- 
tation caused by the accretion of gas from the environment 
of the radio galaxy or by a merger. 

The extrapolation of the radio flux density at 8.3 GHz 
(Rottgering et al. 1994; van Ojik et al. 1996) provides 
expected synchrotron flux densities of 0.2 and 0.4 mJy, 
respectively, at 800 yum and 1.3 mm. We observed this 
source, both with the IRAM telescope and the JCMT, 
but we obtained only non-significant detections at the 2 — 
3-cT level (see Table 1). The IRAM observations were per- 
formed on two different nights. Although the combined 
observations of the two nights provide a formal 3-cr de- 
tection, we found that this result is not reliable because 
the source was detected only during the first night. In fact, 
deeper observations with the IRAM telescope failed to de- 
tect the galaxy and provided a 3-cr upper limit < 1.5 mJy 
(R. Chini, private communication). 

However, we can see that our 3-cr upper limits 
(<9.3mJy at 800 /im; <2.6mJy at 1.3 mm) provide a rel- 
evant result. In fact, the inferred total dust masses are 

< 1.0 x 10 8 and < 1.3 x 10 8 M using the 800-^m and 1.3- 
mm upper limits, respectively. The most stringent limit on 
the dust mass (provided by the JCMT observation) im- 
plies that the amount of dust in 1243+036 is lower than 
that inferred for those high-z radio galaxies and quasars 
so far detected (see Hughes et al. 1997 and references 
therein). The upper limit on can be lowered still fur- 
ther if we use the limit provided by Chini at 1.3 mm, which 
gives Md < 7.7 x 10 7 M Q . It is also important to no- 
tice that CO observations of this galaxy have provided 
a stringent limit on the amount of molecular hydrogen 

< 5 x 10 10 M Q (van Ojik 1995). 

PKS 1251-407 

To date, this is the furthest known radio-loud quasar 
(Shaver et al. 1996). Our SEST observation provides a 
tentative (2.6-cr) detection. Figure 2 shows the SED of 
this quasar. The 1.3-mm upper limit hints that the syn- 
chrotron spectrum steepens at high frequencies, as do 
the other two radio-loud quasars, MRC 1043—291 and 
PKS 1354-107. 

PKS 1354- 107 



The optical counterpart of this radio source was iden- 
tified by Shaver et al. (P.A. Shaver, J.V. Wall, K.I. Keller- 
mann, OA. Jackson, M.R.S. Hawkins, private communi- 
cation) with a quasar at z = 3.0. The available radio flux 
densities (see Figure 2) suggest the presence of a very flat 
spectrum from 2.7 to 8.4 GHz. Our SEST observation pro- 
vides a 3-cr upper limit of <15.3mJy, implying a sharp 
steepening of the synchrotron spectrum at high frequen- 
cies. 

4. The case of MG 1019+0535: a dusty radio 
galaxy 

Although the JCMT observation provides only a marginal 
detection, its combination with the IRAM detection of an 
excess over the synchrotron spectrum strongly suggests 
the presence of thermal dust emission from this galaxy 
(Fig. 3). 

Hughes et al. (1997) have demonstrated that the uncer- 
tainties in calculating the mass of dust responsible for the 
optically thin, thermal, submillimctre emission are: (i) our 
limited knowledge of the rest-frame mass absorption co- 
efficient, Kd, and how this quantity varies with frequency; 
(ii) the dust temperature (7d), and, finally, (iii) the un- 
known values of H and qo. Unfortunately, these problems 
are coupled. For example, Hughes et al. (1993) noted that 
there is a trade off between Td and the critical frequency 
at which the dust becomes optically thick, v . 

Typically, each of these uncertainties can account for 
changes of up to a factor <~ 5 in derived dust masses (see 
Hughes et al. 1997 for more details). However, a conserva- 
tive value for the dust mass can be obtained if the para- 
maters we use to estimate the dust mass are taken such 
that the dust mass is minimised. We further note that 
(since for high-z objects the submillimetre observations 
sample the Rayleigh- Jeans region) the slope of the dust 
spectrum is not a function of temperature. 

Our measurements of MG 1019+0535 at 240 and 
384 GHz suggest that the submillimetre spectral index (a, 
where S v oc v a ) is large and positive (a = +4.2 ± 1.2). 
We recall that the maximum allowed spectral index for 
self-absorption is +2.5. Our result therefore rules out the 
possibility that the emission is due to self-absorbed syn- 
chrotron radiation (Chini et al. 1989). However, given the 
uncertainty of the 384-GHz flux density, data at more fre- 
quencies are needed to better constrain a. 

Strong support for the thermal nature of the submil- 
limetre emission is provided by our deep measurements 
at 22 and 43 GHz using the VLA. These show that the 
steepening centimetre spectral index (Figs. 1 and 3) be- 
comes still more negative as it approaches the millime- 
tre domaine; the predicted contribution at 240 GHz from 
the dominant centimetre component lies several orders of 
magnitude below the measured 240-GHz flux density. 

At first sight this indicates that the frequency depen- 
dence of the dust grain emissivity (or the emissivity index, 
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f3) is +2.2 ± 1.2, which encompasses the range normally 
quoted for interstellar grains (1.0 < (3 < 2.0) as well as 
some less physical values (/? > 2.0). However, the redshift 
is high and the rest-frame frequency of the observed 374- 
GHz emission is close to the turnover of the dust spec- 
trum, so we do in fact require a high value of /3 to fit 
both the 240- and 384-GHz data. For (3 = +2.0, we find 
that 35 < T d /K < 180. The lowest temperatures (35 K) 
are found for an optically thin solution; the highest tem- 
perature (180 K) is permitted when we allow the dust to 
become optically thick (say at vq = 1.5 THz or 200 /xm) 
and to be constrained by the IRAS upper limits. 

Although it is clear that our observations do not con- 
strain stringently the dust temperature, we can make 
use of the usual theory (Hughes et al. 1997) to esti- 
mate the mass of dust responsible for the emission de- 
tected by IRAM and JCMT. For 35 < T d /K < 180, 
0.17 < M d /10 8 M Q < 1.8. 

It is reassuring that T d = 40 K is viable since this is the 
temperature of the dust measured in the z — 4.25 radio 
galaxy, 8C 1435+635 (Ivison et al. 1998). For T d = 40K, 
adopting the same dust parameters as Ivison et al., we 
derive M<j = 2 x 10 8 /igQ 2 M which, when compared with 
the dust mass estimate of 4 x 10 8 /i^ 2 M for 8C 1435+635, 
suggests that powerful radio galaxies evolve significantly 
between z — 4.25 and 2.76 (though we note that the rest- 
frame 6-cm luminosity of 8C 1435+635 is around 5 times 
that of MG 1019+0535 and that observations of complete 
samples of radio galaxies spanning a range of redshifts and 
radio luminosities will be required to trace their evolution 
in detail). 

4.1. Modelling the UV-to-FIR SED 

The interpretation of the observed spectral energy dis- 
tributions (SEDs), is not straightforward, since a non- 
thermal contribution cannot be neglected and the com- 
monly used population synthesis models do not al- 
low for dust extinction. Using the same approach fol- 
lowed by Mazzei & De Zotti (1996), based on chemo- 
photomctric population synthesis models incorporating 
extinction and re-emission by dust and accounting for 
non-thermal emission, we have attempted to analyse the 
SED of MG 1019+0535. We recall here that Dey et al. 
(1995) estimate an upper limit to the internal reddening, 
E B -v < 0.43 mag. 

Our IRAM (and, marginally, JCMT) observations 
strongly favour the presence of dust. Although our spec- 
tral coverage is rather poor, we attempt to fit the overall 
SED of this galaxy, from the optical to 1.25 mm in the 
observed frame, with the aim of constraining the evolu- 
tionary properties of MG 1019+0535. For this study, we 
assume that the source of the submillimetre continuum 
radiation is component A (as suggested by its depleted 
Lya emission) , and we adopt its optical and near-IR fluxes 
accordingly (Dey et al. 1995). 



We have computed several models with Salpctcr's ini- 
tial mass function (IMF) and different lower mass limits, 
mi, as described in Mazzei & De Zotti (1996) (and refer- 
ences therein) . For a given model we derive the age of the 
system which matches the data, with different amounts of 
non-thermal AGN emission. We find an interesting result: 
these data are well matched by models which always cor- 
respond to a 0.8 — 1-Gyr-old host galaxy, accounting for a 
non-thermal contribution ranging from 50 to 90% of the 
total flux density at 0.6 /jm (see Fig. 4). According to this 
result, MG 1019+0535 cannot be considered a "primaeval" 
galaxy candidate because the bulk of its stellar population 
is significantly evolved. For H = 50kms _1 Mpc -1 , the 
formation redshift, Zf mm , of MG 1019+0535 is between 10 
and 4 if q = 0.5 and < 5 if qo = 0. In the following we 
will refer to the first value of qo- 

The expected bolometric luminosity is always larger 
than 10 13 i Q ; in particular, this rises by a factor 2.5 if 
to; = 0.5 Mq. This corresponds to a residual gas fraction 
of 2%, i.e. to a total gas mass of about 4 x 10 10 M Q , which 
is well below the Evans et al. (1996) upper limit, with 
Eb-v = 0.13, a total barionic mass, Mbarion, of around 
1.7 x 10 12 M Q , and a star-formation rate of 800M Q yr _1 . 
In this scenario, the hot stars are almost completely ob- 
scured by dust which, heated by their radiation field, 
transfers their bolometric luminosity to the far-IR wave- 
length regime. Models with lower m; require larger Mbarion 
and higher star-formation rates. We derive Eb-v < 0.26 
for a residual gas fraction as large as 30% — the largest 
allowed by models — and M ba rion < 4.5 x 10 12 M Q . 

The available data can be fully accounted for by 
opaque models like those already used by Mazzei & De 
Zotti (1994) to fit the spectrum of the ultraluminous 
galaxy IRAS F10214 + 4724. However, there is still con- 
siderable latitude for modelling. Crucial constraints may 
be provided by ground-based submillimetre measurements 
and by observations with the Infrared Space Observatory 
(ISO)] these measurements will help to define the shape 
of the far-IR SED, so settling the dust temperature, and 
the role of PAHs in the near-IR spectral range. 

5. Results on radio- loud quasars 

The three quasars show a sharp steepening of their syn- 
chrotron spectra between rest-frame radio and submillime- 
tre wavelengths. It is important to recall that variability 
may affect the flux densities of radio-loud quasars obtained 
at different epochs. In order to estimate the magnitude of 
the spectral break in the millimetre region, we evaluate 
the expected synchrotron flux density at ^ O b s =230 GHz, 
deriving the spectral index from the available data and 
extrapolating the 5-GHz radio flux densities to 1.3 mm. 
If two flux densities are available at the same frequency, 
we adopt the average value. For PKS 1251—407, we use 
only the quasi-simultaneous 1.4 — 15-GHz data of Shaver 
et al. (1996). For MRC 1043-291, we find that a(0.408 - 
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5 GHz) = —0.32, and that the expected flux density 
at 230 GHz is 5 exp (230) - 152 mJy. For PKS 1251-407 
and PKS 1354-107, we find a(1.4 - 15 GHz) = -0.37, 
a(2.7 - 8.4 GHz) = -0.07, with S eX p(230) - 50mJy and 
S'exp(230) ~ 210 mJy, respectively. If we compare the ex- 
pected flux densities with the observed upper limits, we 
find that 5 , cxp (230)/5 , obs (230) > 5 - 14, i.e. the observed 
flux density can be more than an order of magnitude 
lower than expected. We cannot know what fraction of 
this cut-off may be due to variability; however, the typi- 
cal maximum amplitude of the variability of flat-spectrum 
radio sources at A ~ 1mm — 11cm, over timescales of 
months — years, is 40 — 60% (Jones et al. 1981 and refer- 
ences therein). It therefore seems unlikely that variability 
alone can explain the sharp cut-offs observed in the three 
quasars in our sample. 

Moreover, high-frequency spectral turnovers have been 
often observed in radio-loud quasars at lower redshifts 
(Chini ct al. 1989; Antonucci, Barvainis & Alloin 1990; 
Klein et al. 1996). In particular, the spectral break oc- 
curs in radio-loud quasars with either flat or steep radio 
spectra, whereas in BL Lac objects the submillimetre flux 
densities lie on an extrapolation of the radio spectrum 
(Knapp & Patten 1991). Our observations show a turnover 
at logf rest (Hz) < 11.8 — 12, but the lack of further spec- 
tral information does not allow us to infer whether the 
steepening of the synchrotron-sub-mm spectra is related 
to thermal emission from dust peaking at the higher fre- 
quencies. It is important to recall that the question of 
dust in high-rcdshift, radio-loud quasars is very relevant 
to unification models (see Baker & Hunstead 1995 for a 
discussion of dust in radio-loud quasars) and to test the 
hypothesis that many quasars may be missed in optical 
surveys (see Webster et al. 1995). Future observations with 
ISO and SCUBA will extend the SED coverage to higher 
frequencies for a large number of objects. 

6. Concluding remarks 

We have presented observations of the rest-frame far- 
infrared continuum of a small sample of radio-loud AGN 
(4 radio galaxies and 3 quasars) with 2 < z < 4.5. 

One of the main findings is the detection of thermal 
continuum emission from the radio galaxy MG 1019+0535 
(z ~ 2.8). This is particularly important since it confirms 
the suggestion that its weak Lya emission is probably due 
to dust extinction. We attempt to estimate the temper- 
ature of the dust in order to derive the dust total mass, 
but the present data data do not allow us to constrain 
stringently the range of possible temperatures. However, 
we estimate that the total dust mass should be of order 
0.2 - 2.0 x 10 8 M Q for temperatures T d ~ 35 - 180 K. The 
overall rest-frame UV-FIR SED can be accounted for by a 
relatively evolved host galaxy (age ~0.8 Gyr) experiencing 
an episode of vigorous star formation. 



The radio galaxy 1243+036 (z ~ 3.6), on the other 
hand, seems to have a lower amount of dust and molec- 
ular gas than those active galaxies that have so far been 
detected at high redshifts, indicating that the properties 
of the ISM in active objects at z > 2 can be rather inho- 
mogeneous. 

Finally, the three flat-spectrum radio-loud quasars 
(2 < z < 4.5) observed at 1.3 mm show evidence of a spec- 
tral turnover at high frequencies. Our data do not allow 
us to understand what fraction of the turnover is due to 
the effects of variability, but we conclude that variability 
cannot be the only cause. 
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Fig. 1. Spectral energy distributions of radio galaxies. The flux densities and the relative references are listed in Table 3. 



Fig. 2. Spectral energy distributions of radiodoud quasars. The flux densities and the relative references are listed in Table 3. 
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Fig. 3. Spectral energy distribution of MG 1019+0535. The lines drawn in the rest-frame far-IR are modified 35- and 180-K 
blackbodies, with a (3=+2 frequency dependence for the dust-grain emissivity, representing the most extreme dust temperatures 
compatible with our data. The 180-K blackbody is optically thick at 200 /xm. Key: circles — VLA, JCMT, IRAM and IRAS 
measurements described in the text; diamonds — measurements from Dey et al. (1995); squares — measurements from Griffith 
et al. (1995), Becker, White & Edwards (1991), Wright & Otrupcek (1990) and White & Becker (1992). 
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Fig. 4. The fits to the SED of MG 1019+0535 as obtained with the Mazzei & De Zotti (1996) model. Thick lines correspond 
to the overall match of the SED of MG 1019+0535 for models including a non thermal contribution, stellar emission and dust 
effects (see text); a) shows the results for two models corresponding to a host galaxy 1-Gyr old, short-dashed (m; = 0.01 Mq) 
and long-dashed lines (mi = 0.5 Mq), with a non-thermal contribution at 0.64 /im of 70 and 80% respectively (thin line); the 
overall match for a host galaxy 0.8-Gyr old (dot-dashed line, mi — 0.1 Mq) with a non-thermal contribution of 50% at the same 
wavelength is also shown; in b) are the results for the same models raising the non-thermal contribution at the same wavelength 
to 90%; short-dashed curves require a dust temperature of 46 K instead of 60 K. 
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